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Abstract 

In the two-dimensional isotropic parabolic potential barrier V(x,y) = Vq — 
m'y 2 (x 2 + y 2 )/2, though it is a model of an unstable system in quantum mechanics, 
we can obtain the stationary states corresponding to the real energy eigenvalue Vo- 
Further, they are infinitely degenerate. For the first few eigenstates, we will find 
the stationary flows round a right angle that are expressed by the complex velocity 
potentials W = ±jz 2 /2. 



1 Introduction 



It is well known that the two-dimensional harmonic oscillator is equivalent to the dynam- 
ical system consisting of the two independent one- dimensional harmonic oscillators — the 
energy eigenvalues of the two-dimensional oscillator are given by the sum of the energy 
eigenvalues of the one-dimensional oscillators and the eigenstates of the system are given 
by the product of the eigenstates of the one-dimensional oscillators. When degenerate 
eigenstates of the two-dimensional oscillator are superposed with a suitable weights, the 
new states will be the eigenstates of orbital angular momentum. These results were stud- 
ied a long time ago by Dirac Jl||. 

In the present paper we will investigate the two-dimensional parabolic potential bar- 
rier, which is a model of an unstable system in quantum mechanics, on the same lines 
as the two-dimensional harmonic oscillator. This model is equivalent to the dynamical 
system consisting of the two independent one-dimensional parabolic potential barriers. 
The one-dimensional potential barrier was studied by the authors |2|, ||| . It is shown that 
the energy eigenvalues are complex numbers and the corresponding eigenfunctions are 
expressible in terms of the generalized functions of a Gel 'f and triplet. In two dimensions 
the exact solutions of the eigenvalue problem of this model separate into four types. We 
will take the types of two of the four types in § [3.1| , for which the solutions are expressed 
by generalized eigenfunctions belonging to complex energy eigenvalues and represented 
by diverging and converging flows. In these two types the solutions will also be the eigen- 
states of orbital angular momentum. We will study the other two of the four types in 
§ p.2| . In these two types all the solutions are infinitely degenerate and involve the special 
solutions with real energy eigenvalue. Such special solutions are represented by stationary 
corner flows. From the hydrodynamical point of view they can be described by some kind 
of quantum "velocity" and complex velocity potential discussed in a previous paper f|] . It 
has been pointed out to us by R. Jackiw that this velocity was originally introduced half 
a century ago by Madelung ||. Such a velocity is still useful in present-day high-energy 
physics 0. We shall see that, for the first few solutions, the velocities reviewed in § g are 
solenoidal, so the corresponding complex velocity potentials must exist. These complex 
velocity potentials for the two-dimensional parabolic potential barrier describe the flows 
round a right angle. 

2 Complex velocity potentials in quantum mechanics 

In the present section we shall summarize the features of the velocity in quantum me- 
chanics. 

Define the velocity of a state ip(t, r) in non-relativistic quantum mechanics by 
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where j(t,r) is the probability current 



j(t,r) = ®[il>{t,r)*{-ihV)if)(t,r)]/m, (2.2) 

and m is the mass of the particle. 

If all variables are separable, this velocity is in general irrotational, namely, the vor- 
ticity defined by 

u> = V x v (2.3) 

vanishes. Then the velocity in irrotational flow may be described by the gradient of the 
velocity potential <P, 

v = W. (2.4) 
We now proceed to study only the two-dimensional flow. Let us consider the velocity 



(2TTJ) which is solenoidal, namely 

V .„^ + ^ = 0. (2.5) 
ox ay 

The velocity in two-dimensional flow satisfying this solenoidal condition may be described 
by the rotation of the stream function 

v * = w Vy = ~ ^ ■ ( 6) 



Further, in the two-dimensional irrotational flow, equations ( |2.4| ) and ( |2.6| ) can be 
combined into Cauchy-Riemann's equations between the velocity potential and the stream 
function. We can therefore take the complex velocity potential 

W(z)=$(x,y)+i&{x,y), (2.7) 

which is a regular function of the complex variable z = x + iy. For example, the flow 
round the angle 7r/a is expressed by 

W = Az a , (2.8) 

A being a number. With a — 1, this expresses the uniform flow of the two-dimensional 
plane wave in quantum mechanics. There are some elementary examples of complex 
velocity potentials which express the two-dimensional flows in quantum mechanics ||. 

The above-mentioned method will be applied to the stationary flows of the two- 
dimensional parabolic potential barrier in 
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3 The parabolic potential barrier in two dimensions 

The Hamiltonian of the two-dimensional isotropic parabolic potential barrier is 

6 =-H& + w) +Vi -7*V + rt- (3-D 

where Vq G R is the maximum potential energy, m > is the mass and 7 > is propor- 
tional to the square root of the curvature at (x, y) = (0, 0). 

A state is represented by a wave function U(x, y) satisfying the Schrodinger equation, 
which now reads, with H given by (|3.1| ), 

( + |l) "fa V) + {V0- \W (x 2 + y 2 ) } U(x, y) = EU(x, y). (3.2) 

The energy eigenvalues of (|3.2|) will be the sum of the energy eigenvalues of the one- 
dimensional parabolic potential barrier in the ^-direction and ^-direction, respectively, 
i.e. 

E nxny = E nx + E Hy (3.3) 

and the eigenf unctions belonging to these energy eigenvalues will be the product of their 
corresponding eigenf unctions 

Un x n y (x,y) =u nx (x)u ny (y). (3.4) 

With the notation of preceding papers 0, |3|, the energy eigenvalues of the one- 
dimensional parabolic potential barrier are 

< = ^o^(n 9 + 0fry K = 0,l,2,---) (3.5) 

and the corresponding eigenf unctions are 

<(?) = e ±l ^ 2/2 H±((3q) (0 ee ^Jh) , (3.6) 

where H^ {f3q) are the polynomials of degree n q , and the numerical coefficients are 
discarded. The eigenfunctions are generalized functions in 1 S(1R) X of the following 
Gel'fand triplet, 

S(R) c L 2 (R) c S(R) X , (3.7) 

where L 2 (M) is a Lebesgue space and S(R) is a Schwartz space. The work also shows 
that the index + means only outward moving particles and the index — means only inward 
moving particles. 
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Thus the results ( |3.3| ) and ( |3.4| ) of the two-dimensional parabolic potential barrier 
separate into four types: 

Type 1. E++ y = E+ x + E+ =V Q - i(n x + n y + 1)%, 

UtX&y) = <»<» = e + ^ {x2+y2)/2 Ht([3x)H+ y ((3y). 

Type 2. E+~ ny = + E~ y = V - i(n x - n y )h 1} 

Ut- y (x,y) = ul(x)u- y (y) = e^^H^H^y). 

Type 3. E~+ ny = £T + E+ y = V + i{n x - n y )H-y, 

u~X^y) = <»<» = e-^ 2{x2 ~ y2)/2 H- x ((3x)H^ y (p y ). 

Type 4. E~ y = E~ x + E~ y = V + i{n x + n y + 1)%, 

U-- y (^y) = u- n Sx)u- y {y) = e-^ 2+ y^H- x ((3x)H- y (f3y). 

These eigenfunctions > Un~n y -> ^n^n y are a l so generalized functions in iS(IR 2 ) x 

of the Gel'fand triplet 

S(R 2 ) C L 2 (R 2 ) c 5(R 2 ) X (3.8) 

instead of ( |3.7| ). Note that the eigenfunctions of types 4 and 3 are conjugate complex 
functions of types 1 and 2, respectively, i.e. 

and 

U%,{x,vY = UZ%(x,y). 
3.1 Diverging and converging flows 

Let us consider first the types 1 and 4. In this case the energy eigenvalues E^ y are 
always complex numbers and the time factors corresponding to them are 

e -iE^„ y t/h _ e -iV t/h e ^(n x +n y +l)-yt 

Thus the solutions of type 1 are well defined when t > 0, and those of type 4 are well 
defined when t < 0, according to the time boundary condition that time factors of an 
unstable system are square integrable f2|. Also, U^ y (x, y) represent particles moving 
outward from the center as in fig. [I], and (x,y) represent particles moving inward 
to the center as in fig. ||. Thus we shall call these types diverging and converging flows, 
respectively. Note that a time reversal occurs resulting in the interchange of the diverging 
and converging flows. 
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FlG. 2: Converging flows. 



For n x = n y = 0, we get the energy eigenvalue 

E±± = V Tih 1 (3.9) 

and only one eigenfunction 

O^) = ^ 2(x2+y2)/2 , (3-10) 
respectively. For n x + n y = I, namely n x = 1, n y = and n x = 0, n y = 1, we get 

^ = Eo? = VoT 2ihry (3.11) 



(3.12) 



and two eigenf unctions 

U^ ± (x,y) = 2f3xe ± ^ 2+ ^ 2 , 
U± ± (x,y) = 2Pye ±l ^ x2+ y 2 V 2 . 

There is twofold degenerate state of types 1 and 4 with n x + n y = 1. For n x + n y = 2, 
namely n x = 2,n y = 0; n x = 1, n y = 1; ra x = 0, n y = 2, we get 

^ = ^ = = V T 3ifry (3.13) 

and three eigenfunctions 

C^(x J y) = (4^x 2 =F2i)^* a+ * a >/ 2 I 

t/±fc( x , y ) = Af3 2 xye ±l ^ x2+ y^ 2 , } (3.14) 
U^y) = {^ 2 y 2 T2i)e^ 2+ y^ 2 . 

There is threefold degenerate state of types 1 and 4 with n x + n y = 2. Generally, there is 
in + l)-fold degenerate state of types 1 and 4 with n x + n y = n. This result is just the 
same degree of degeneracy as the two-dimensional harmonic oscillator. 

For the further discussion of the state of types 1 and 4, we now pass from the Cartesian 
coordinates x, y to the two-dimensional polar coordinates r, ip by means of the equations 



x = r cosy?, 
y = r sin (p. 



(3.15) 



If in the new coordinates we superpose above-mentioned eigenstates with suitable weights, 
the result will be the eigenstates of orbital angular momentum L defined by 

L = -iH [x%- - y4-] = -if*S-- ( 3 - 16 ) 
\ ay ox J dip 

For n x = n y = 0, the eigenfunction Q3.10Q will be 

U±±(r,p)=e^ 2 / 2 . (3.17) 
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Thus Uqq (r, ip) is independent of if and has zero orbital angular momentum. For n x +n y = 
1, a linear combination of the eigenfunctions Q3.12j ) give 

U±±(r, <p) + iU^(r, <p) = 2(3re^ 2 /^^ 1 

These states are eigenstates of L with eigenvalues H and —H, respectively. For n x + n y = 2, 
(gUD give 

^(r, + 2iE/£ t (r, ¥») - ^(r, *>) = 1 

^(r, ¥>) + t^(r, *>) = 4(/3 2 r 2 =f Oe^' 2 , I (3.19) 
U±±(r, <p) - 2iU?±(r, tp) - U±±(r, <p) = ^e^l\~^. J 

These states are also eigenstates of L with eigenvalues 2H, 0, and —2h. There is a similar 
procedure for large values of n x + n y . 

3.2 Corner flows 

Let us now study the types 2 and 3. In this case the energy eigenvalues are also 

in general complex numbers, but with the striking difference that all the eigenstates 
belonging to each energy eigenvalue are infinitely degenerate. The corresponding time 
factors are 

Thus the solutions of type 2 are well defined when t > 0, and those of type 3 are well 
defined when t < 0, for the case of n x > n y , and vice versa. Also, U^~ ny {x,y) represent 
particles which, coming from the y-direction, round the center and, go off to the x-direction 
as in fig. ^ and U~^ y (x, y) represent particles which, coming from the ^-direction, round 
the center and, go off to the y-direction as in fig. [|. Note that a time reversal occurs 
resulting in the interchange of these corner flows. 

Stationary flows 

The above time factors now show that for the case of n x = n y , there are stationary 
flows. For = 0, 1, 2, • • • , the energy eigenvalues associated with stationary 

flows are the same real number: 

E% = V . (3.20) 
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FlG. 3: Corner flows moving from the y-direction to the re-direction. 




FlG. 4: Corner flows moving from the x-direction to the y-direction. 
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The first few infinitely degenerate eigenf unctions belonging to this energy eigenvalue ( |3.20| ) 
are 



.,2 



U^(x,y)=4f3 2 xye ± ^ x -y ^ 2 , 
U?*(x, y) = 4 [4f3*xV + 1 ± 2^ V - y 2 )} e^~^\ 



(3.21) 



For the further study of the stationary flows of the two-dimensional parabolic potential 



barrier with the Hamiltonian (|3.1|) , it is convenient to make a transformation to the 
rectangular hyperbolic coordinates u,v, given by 



2 2 

u = x — y , 



The eigenf unctions 



v = 2xy. 

will become in the new representation 



(3.22) 



u,v) 



±i/3 2 u/2 



U^(u,v) = 2f3 2 ve ±ip2u/2 , 

U±*(u, v)=A (/?V + 1 ± 2if3 2 u) e^ 2u ' 2 , 



(3.23) 



The factors e ±l/3 u l 2 occurred in Q3.23| ) describe plane waves in the w-plane, i.e. the 
motion of the wave e t/3 U//2 is in the direction specified by fig. [3] and that of the wave 
e -i/3 u/2 j g j n direction specified by fig. [|. These eigenf unctions substituted in fl2.2p 
give the rectangular hyperbolic coordinates j^nm 3nnv °f 3nn-> which are the probability 
currents of the states U^. They give 



;±=f. 



u,v) 



±lh u /2, j^(u,v) = 0, 

±2 1 (3 4 v 2 h u , jt£(u,v)=0, 

±87 {(/?V + 5)(/?V + 1) + 4/3 V} h u , 

^QA^f3 A uvh v , 



where the scale factors h u = h v = 2\/u 2 + v 2 . Thus can never depend on the time t. 
We see from this result the suitability of the term "stationary flows" . 



The flows round a right angle 

The above probability currents now show that for the case of n = and 1, there are 
stationary flows which move along the hyperbolas (each line with v constant). To get an 
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understanding of the physical features of this flows it is better to work with the velocity 
defined by (|2.1| ). For n = and 1, the velocities give the same result 



v^ = ±- 1 h u , v^ = 0. 

Taking the rotation of them, the vorticity ( |2.3|) becomes 

d ( t> ±=F \ d ( v^ 



UJ 



du \ h 



dv V h 



0. 



These equations will hold generally in quantum mechanics M, and therefore the velocity 
potentials defined by ( |2.4| ) must exist. If we transform to rectangular hyperbolic coordi- 
nates u, v, equations (|2.4|) become 



h d<P 

V u "u"^, j V v 

du 



h 



d<P 

' dv ' 



(3.24) 



and the velocity potentials for n = and 1 are thus 



2 ' 



(3.25) 



Note that they are proportional to the phase factors of (|3.23|) . Further, the divergence 

(ED g ives 

d ( v , ±=F \ d ( v , ±=F ~ 



du \ h 



dv \ h. 



0. 



Thus the velocities are solenoidal for n = and 1, so that we can obtain the stream 
functions defined by (|2.6|). The equations (|2.(j|) are also expressed, as in equations (^.24|), 



(3.26) 



h — h — 

OV OU 

and the stream functions for n = and 1 are thus 

2 ' 



(3.27) 



For the states represented by the first and second of equations ( 3.21 ) or (|3.23|) , the complex 
velocity potential ( |2.7| ) gives, from ( ^.25| ) and ( |3.27| ) 



1 i 
W + = ±-7« ± -jv 

1 9 



(3.28) 



since z 2 = u + iv. Equations ( |3.28|) are of the form (|2.8| ) with a = 2, and they show that, 
forn = and 1, the complex velocity potentials of the two-dimensional parabolic potential 
barrier express the flows round a right angle. One could work out in terms of Cartesian 
coordinates and one would be led to the same conclusion. 
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4 Discussion 



We have obtained the exact solutions of the two-dimensional parabolic potential barrier. 
One class of the solutions is diverging and converging flows of § |3.1| . These solutions 
are always complex energy eigenvalues and generalized eigenfunctions, which mean that 
the diverging and converging flows are not stationary. These generalized eigenfunctions 
can be obtained from the eigenfunctions of the two-dimensional harmonic oscillator by 
the analytical continuation, in the same way as the one-dimensional parabolic potential 
barrier 0. Again, they can be superposed to give the eigenstates of orbital angular 
momentum. For these solutions, however, the method mentioned in § [2| was not applicable. 
As an example we try to calculate the divergence ( |2.5|) for the states ( |3.10| ) or ( |3.17| ). The 
result is 

V • vff = ±2 7 ^ 0. 

Now the ± sign shows that Uq + is connected with the diverging flow and Uq ~ is connected 
with the converging one. Thus the velocities of diverging and converging flows cannot be 
solenoidal. Therefore the stream functions and the complex velocity potentials do not 
exist. This result still holds for large values of n x + n y . 



The other class of the solutions is corner flows of § |3.2| . All the solutions are infinitely 
degenerate and involve the stationary flows with real energy eigenvalue. It should be noted 
that there are no stationary flows in the one-dimensional or three-dimensional isotropic 
parabolic potential barrier. For n = and 1 in the stationary flows, we have found the 
flows round a right angle that are expressed by the complex velocity potentials ( |3.28|) . 



But for n ^ 2, the complex velocity potentials do not exist, because the imaginary parts 
of the polynomials H^(/3q) in (|3.6| ) cause the stream lines to depart from hyperbolas. 
One may, however, find a new flow as the result of a kind of superposition of the infinitely 
degenerate states. 

One would expect to be able to get a more direct solution of the eigenvalue problem 
of the Hamiltonian ( |3.1| ) by working all the time in the two-dimensional polar coordi- 
nates, instead of working in the Cartesian coordinates and transforming at the end to the 



two-dimensional polar coordinates, as was done in § |3.1| . But under suitable boundary 
conditions in the two-dimensional polar coordinates, one would obtain only diverging and 
converging flows of § |3TT| , i.e. the lack of corner flows of § |37^. It is also pointed out that 
one can get the only diverging and converging flows from the analytical continuation of 
the solutions of the two-dimensional harmonic oscillator. These facts mean that the choice 
of coordinate systems is quite important in the eigenvalue problem of the unstable system 
in non-relativistic quantum mechanics, since coordinate systems impose a restriction on 
the symmetry of boundary conditions. The source of the conclusion lies in the existence 
of a very large class of solutions for the unstable system. 

It is rather surprising that such a stable idea as stationary flows should appear in the 
parabolic potential barrier in this way. Actually we have shown in another paper |7j that 
the dynamical system composed of several of these potential barriers forms the quasi- 
stable semiclassical system. 



11 



References 

P. A. M. Dirac, The Principles of Quantum Mechanics (Clarendon Press, 1930). 

T. Shimbori and T. Kobayashi, Complex Eigenvalues of the Parabolic Potential Barrier 
and Gel'fand Triplet, math-ph / 99 10009] , to appear in Nuovo Cim. B. 

T. Shimbori, Operator Methods of the Parabolic Potential Barrier, |quant-ph / 9912075 . 



T. Shimbori and T. Kobayashi, ''Velocities" in Quantum Mechanics, |quant- 
ph/0004086 . 



E. Madelung, Z. Phys. 40 (1926) 322. 

For example, D. Bazeia and R. Jackiw, Ann. Phys. 270 (1998) 246. 

T. Kobayashi and T. Shimbori, Statistical Mechanics for States with Complex Eigen- 
values and Quasi-stable Semiclassical Systems, |cond-mat / 0005237 . 



12 



